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Abstract 

NMR imaging without any radio-frequency magnetic field is explained 
by a quantum treatment of independent spin The total magneti¬ 
zation is determined by means of their individual wave function. The 
theoretical treatment, based on fundamental axioms of quantum me¬ 
chanics and solving explicitly the Schrodinger equation with the kinetic 
energy part which gives the motion of free spins, is recalled. It explains 
the phase shift of the spin noise spectrum with its amplitude compared 
to the conventional NMR spectrum. Moreover it explains also the rela¬ 
tively good signal to noise ratio of NMR images obtained without a RF 
pulse. This derivation should be helpful for new magnetic resonance 
imaging sequences or for developing quantum computing by NMR. 
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Introduction 


The work of Norbert Muller and Alexej Jerschow [lj on nuclear spin noise 
imaging has reactivated the fundamental idea given by Felix Bloch j2] in 1946 
that N nuclei with magnetic moment /j can be measured because of “statis¬ 
tically incomplete cancellation”. A recent review on nuclear spin noise which 
corresponds to this idea can be found in reference j3|. The aim of this work 
is to set the frame for a quantum mechanical description of an ensemble of 
N identical non interacting spins \ at finite temperature which experience the 
presence of main static and radio frequency (RF) magnetic fields in magnetic 
resonance experiments. We will recall the elementary quantum description 
of such a system which explains the small induced current by an ensemble of 
spins at equilibrium and the possibility to make NMR images without any 
RF magnetic field only with time dependent gradient magnetic fields. This 
derivation should be helpful for new magnetic resonance imaging sequences 
or for developing quantum computing by NMR. 

1 Theoretical description 

1.1 One spin | in a static magnetic field 

For one nuclear spin \ in a B 0 magnetic field, the potential part of the 
Hamiltonian is: 


Jrf? — ft ■ Bq — — 7 / • Bq, 


( 1 ) 


where 7 is the gyromagnetic ratio of the nuclear spin and / is the spin that 
verifies the following equation: 



( 2 ) 


with <t a vector defined by the following Pauli matrices: 



( 3 ) 


Matrix notation of equation (|T|) is given by: 
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By convention, the vector B 0 = (B x , By, Bz) is placed in the Oz-direction, 
that is B 0 = (0, 0, B 0 ) and the xOy-plane is called the transverse plane. 
Therefore the matrix form of the Hamiltonian reduces to: 





(5) 


The two eigenvalues of the Hamiltonian Jt?, which give the energy of the 
quantum states, are: E_ = — hBo and E + = ly/hBo which have the fol¬ 
lowing two corresponding eigenvectors: 


-> 




( 6 ) 


The wave function of one spin| at time t — 0 is: 


^(0)) = xi |-) + x 2 |+) = rqe^ 1 |-) + r 2 e lct>2 |+) 


/ r2 g*02\ 


with r\ + T'l 


1 


(7) 


The 1 indice corresponds to the low energy or ground level. Time evolution of 
this quantum system (far from the speed of light) is given by the Schrodinger 
equation: 

MU .*|,( t )>. (8) 

If we write 

h (t) > = (*$]) • < 9 > 

then we have the following differential equations: 


ih 


(x 2 {t)\ 



1 k ( UQX 2 {t)\ 

(t)J 


fx 2 (t)\ 


( 10 ) 


where 


uj 0 = -7H0 


( 11 ) 
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is the resonance or Larmor frequency. The following time dependent wave 
function is the obvious solution of this system of equations: 




\^ne +i (^ +01 ) J 




( 12 ) 


= £((n 0 ,t)|^(0)) 


which dehnes the spin evolution time operator E(ojQ,t). 


1.2 Effect of a RF magnetic field 

In NMR, a transition between the two states |—) and |+) is obtained by a 
B i magnetic held rotating in the transverse plane. B\ is assumed to lie along 
the Ox-direction. This held is generated by an electromagnetic RF wave ob¬ 
tained by an oscillating current in a solenoid surrounding the spin sytem. In 

this case the Hamiltonian has non-diagonal elements due to the RF magnetic 

—* —* 

held B i rotating around B o with an angular velocity lu: 

(13) 

In this case time evolution is no longer trivial. We have the following differ¬ 
ential equations: 

ihd\^ F (t)) ( x 2 (t)\ 

_i h ( u 0 wie l0jt \ _ (x 2 (t)\ 

2 \uie iut -u 0 J [x^t))’ 

where 


LUi = —'fBi. 


(15) 


In the time-independent case (cai = 0) these equations reduce to equa¬ 
tions (10). For the time-dependent case we need to solve the following dif¬ 
ferential equations: 


Ar 2 (tA _ /Wr 2 (t) + u i e luJt xi(t)\ 

2 \uj l e lU}t x 2 {t) -ujQXiit) J ' 
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To solve these two equations we make the following substitutions: 


(17a) 

(17b) 


.wo, 

p(t) = x 2 (t)e l 2 , 
q(t) = X! 

Therefore we have: 

p(t) = (^~ l - {uox 2 (t) +(Jie~ lut x 1 (t)) + l ^-x 2 {t) S j e 1 ^ (18a) 

= - l ^-e~ iuit xi (t)e*^ = --^g(f)e i(a;o_aj)t , 

and 

q(t) = l - (uie lut x 2 (t) -cu 0 Xi(f)) - ^aq(t)^ (18b) 

= = -^p(*)e“ i(w °“ w)t . 

When we take the second derivative of p(t) we obtain the following second 
order differential equation: 

p = - + i(cu 0 - ta)p 

= + z(u 0 - u)p (19) 

= i(u 0 -uj)p- ^p. 

Let A ± be the solutions of the equation 

L0 2 

\ 2 -i(u 0 -uj)\ += 0 . ( 20 ) 

We have 

A± = i ——— with fl = uq — u A = \Jvl 2 + (jj\ (21) 

and 

pit) = C ie A +* + C 2 e x ~\ (22) 
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therefore: 


icopt 

x 2 (t) = p(t)e~ 2 

= (Cie x+t + C 2 e x ~ t ) e~~^~ (23a) 

.<£. / .a a \ 

= e *2* Cie* 2* + C 2 e ! 2 ( , 


with initial value 


x 2 (0) = r 2 e^ 2 . 

A similar computation for x\(t) leads to: 

iuit 


xi(t)=e 2 ( C^e* 2 4 + C 4 e * 2 4 


xi(0) = rie l(t>1 


(23b) 


(24a) 

(24b) 


As x(t),y(t) verify initial values and the differential equations (16) we have: 

(25) 


Ci + C 2 = r 2 e^, 
(A + fl)Ci + CJ 1 C 3 = 0, 


(A 


C 3 + C 4 = ne** 1 , 

n)c 2 - u> iC 4 = 0 


Computation of the constants C % is now straightforward and gives: 


Cl = § ((1 - g)^ 2 - fr ie ^) , (26a) 

C 2 = §((1 + g)r 2 e^ + f r ie ^) , (26b) 

C 3 = | (-fr 2 e^ + (1 + g) ri e^) , (26c) 

C 4 = | (f r 2 e^ 2 + (1 - g)r ie ^) . (26d) 


These expressions are simplified if the angular velocity u of the RF magnetic 
held verifies the resonance condition u = loq and in this case we have A = oji, 

n = o. 

In general case we can write the evolution of the spin wave function in matrix 
notation as: 


IV’rfW) = A(u>, u 0 , ui u t) IV’i(O)). (27) 
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where A matrix is 


A 


—i 77 1 

e 2 a 


e l 2 % 


e-^b] 
e*2 t a ) 


(28) 


with 


a = cosf-fsinf b = sin f (29) 

where a is the complex conjugate of a. In this equation it is possible to 
separate the evolution due to the RF frequency oj. We then obtain: 


A{uj, ujq, cji, t ) = E(u, t ) ■ R(cu, co 0l coi, t ), (30) 


or 


A = 



0 \ (a b\ 

ef ) \ b V 


At the resonance frequency, the matrix A reduces to: 


Ao(^o, 0J1A) 



cos 


UJ\t 


-i sm 




—I sin 




COS 


U)\t 


(31) 


(32) 


Effect of an RF pulse is usually described as a rotation ||3jJ with an angle of 
9i = 0 Jit in the spin space around vector u = (u x ,u y ,u z ) as: 



( cos ^ — iu z sin y (—iu x — u y ) sin ^ 
(—iu x + u y ) sin y cos y + iu z sin y 


(33) 


If u lies in the transverse plane: u = (1, 0, 0) where ^-direction is given by 
the phase of the RF pulse (in other words by the beginning of the RF pulse) 
and therefore the rotation matrix becomes: 


R 


1 

2 


(1,0,0),0i 


cos 

—i sin %■ 


-i sin y 
cos y 


(34) 


We notice that this is exactly the transformation matrix that we have com¬ 
puted at resonance in equation (32), except for the fact that the time evo¬ 
lution E(oj,t ) is missing in this equation. Therefore in the general case, 
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equation (28) gives the evolution of the wave function and equation (30) 


shows that the evolution of the spin system is a time evolution given by ma¬ 
trix E(oj,t) which depends on the RF frequency followed by a rotation of 
angle 9 1 — At — \J (cuq — w) 2 + uf t around vector: 


'Ujui . Ucr 




) 'U"y 0 ? 


(35) 


u, — 


\J (cuo — cj ) 2 + uj\ 

UJq — OJ 

\J (cuo — to ) 2 + uf 

ft is important to note that this is a lion-commutative product. In other 
words, the spin evolution in presence of an RF held is, according to equations 


(28) and (30), a matrix product of a time evolution operator and a rotation 


operator. The wave function after the pulse is given by: 


| V’rfW} = A-\if> 1 (0)) = E-R- I^i(0)> (36) 

which is the same as the equation derived in Abragam textbook [3] but here 
we have calculated explicitly the A matrix. Therefore at resonance the wave 
function at the end of the RF pulse is equal to the wave function multiplied by 
the rotation matrix of angle Q\ = cot which is then multiplied by the evolution 
matrix but there is no time sequence between these two operations. If this 
elementary description of the spin dynamic retrieve the main characteristic 
of NMR experiments, it is not complete. 


1.3 Kinetic part of the Hamiltonian 

As an Hamiltonian corresponds to the sum of the potential and the kinetic 
energies of the physical system under investigation, the potential energy (— jl- 
Bq) must be completed by the kinetic energy (p 2 /2m) or more precisely 
K' = -^ 1 V 2 (F 1 ('r, t)) where F is a function of space and time dependent on 
the spin system. For our purpose we consider only the time dependence K'(t) 
where K' is a scalar dependent on the absolute temperature for a macroscopic 
sample. It is well known that such a system obeys always the Schroclinger 
equation. The complete Hamiltonian for one spin is: 




K' -ft-B 0 


\h 


, , i 2 K' 

w o + ~r 

2 K' 
h 


2K' 

h 

-UJ 0 + 



(37) 
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and the Schrodinger equation becomes: 

fX2r(t)\ = _i_ ((<*>o + K)xzr(t) + Kx 1 T (t)\ 

\xi T (t)) 2 \Kx2T(t) - (a ; 0 - K)x 1 T (t)) ' 

with K = 2 K'/h. 

Here to solve these equations we make the following substitutions: 


(38) 


.OJQ+K , 


. cjq-K . 

q(t) = Xi{t)e 1 2 . 


p(t) = X 2 (t)e l “2 \ 

which gives (if Uq and K are time independent): 

K 2 K 2 

P(t) - iw«p(<) + -j-p(i) = 0 , q(t) + m 0 <j(f) + —?(«) = 0 


(39) 


(40) 


These equations have the same solutions as equation (20) with hi = luq and 
oji = K. The solution of these equations is: 

A-K. 


(C l e i ^ t + C 2 e 


.A -K 




(41) 


A A \ 

cv 2" + C 2 e~ % 2 * 


• A. .A 

C 3 e*2* + C 4 e *2*y 


where Ci, C 2 , C 3 , C 4 are identical to the coefficients obtained in the previ¬ 
ous case and given by equations (26) with exchanging by oj 0 and A by 
y/cug + K 2 . This solution is again the product of an evolution matrix with 
a rotation matrix: 


\Mt)) = E T -R-\Mu)) 

where Ex is a matrix and not a scalar, and R is given by 


R = 


cos ^ 


sin ^ 


-if sin f 


"if sin f 
cos + i ^ sin ^ 


(42) 


(43) 


Therefore the evolution of the wave function of one spin is given by the 
product of the spin evolution operator and a spin rotation of angle 6 = 
a/cUq + K' 2 t around u(K/ A, 0, u 0 /A). As the evolution operator is a rotation 
there is no oscillation between the two states, and without RF or static 
magnetic time evolution, it seems that no NMR signal can be detected. 
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1.4 General description of one spin | 

In the general case the Hamiltonian of a spin | is dependent on the time 
evolution of the magnetic field B(t) and the kinetic energy therefore 

it can be written as: 




1 f. f lB z + K' 7 B x ~ infB y + K'\ 

2 [yB x + ilB Y + K’ - 7 B Z + K' ) 


and the Schrodinger equation becomes: 

/ _ _i_ ((ujz + K)x 2 + (u x - iwy + K)xi\ 
V-ti ) 2 \(ux + iu Y + K)x 2 - (cn z - K)xiy 


(44) 


(45) 


with K = 2K'/h, (jJx = —'yBx, wy- = — 7 By , = —'fBz ■ 

We look for an adequate exponential substitution p(t) = x 2 (t)e A ^ and g(t) = 
Xi(t)e S( ' t - ) which derivation gives: 


p = p 


A-i(u z + K) 


q = q 


B + i(u z ~ K) 


f (ary ~ iojy + K)qe^ A B) 
|(c Ox + iwy + K)pe( B ~ A) 


(46a) 

(46b) 


For A = 7}(lu z ( 0) + K(0))t, B = 0) — li(0))t these equations simplify 

but the second derivation introduce a coupling between the two states: 


P — — \ (U 2 X + 4 + K 2 + 2 Ku x )p + p(d - -B) (47a) 

- f (dry - fchy + k)qe^ A ~ B) 

q = — \ikx 4* ^y T A 2 + 2Iiujx)q + p(B — 4) (47b) 

- §(dry - ichy + K)pe^ B ~ A) 

The constant solution of these equations corresponds to: 

p - icozP + K^x + 4 + K 2 + 2Kuj x )p = 0 (48a) 

q + iu>zQ + j(^x T o>y + A 2 + 2Ku>x)q = 0 (48b) 


In case of constant magnetic held B(0,0,Bo) we retrieve equations (40) and 
the comparison with these equations shows that the general derivation intro¬ 
duces a coupling of the kinetic energy and the X component of the magnetic 
held in case of non constant magnetic helds. The solutions of equations (48a) 
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is, as usually, obtained by the constant solution then varying the constants. 
The constant solution of these equations, equivalent to equation (20), is : 

x 2 (t) = p{t)e~ m = (C ie x+t + C 2 e x - 1 ) e -^+i<)t ( 4 9a ) 

Xl (t) = q(t)e~ m = (C' 3 e _A_t + C 4 e~ x+t ) e \^ z ~ K)t (49b) 


with 


A± = 
A = 


i ( u z ± A) 


^ 1 , 


(50) 


-t~ ojy T K T 2 K lO x 


u — U) r 


where Ci coefficients are given by equations (26) with Q = ujq 
It is important to notice that c o z is now the frequency difference with the 
reference frequency loq which is the resonant frequency at time t — 0. The 
explicit constant solution at each time t is: 

x 2 (t) = r 2 (cos y£cos(-yf + fa) - y sin yf sin(y£ - fa)) (51a) 

+ sin yt sin(—yf + fa) 

- i (rq^sin ytcos(-yt + fa) + r 2 (^ sin yfcos(yf - fa)) 

- i (cos yf sin(— y f + fa)) 

Xi(t) = rq (cos yf cos(— yf + 0i) — y sin yf sin(— yf + <4)) (51b) 

- r 2 y sin yf sin(yf - fa) 

+ i (^i(y sin yf cos(—yf + 0i) + cos yf sin(-yf + 0i))) 

- i (r 2 y sin ^|f cos(yf - </> 2 )) 

which verify |aq(0)| 2 = r 2 , |x 2 (0)| 2 = r 2 and |aq(f)| 2 + \x 2 (t)\ 2 = 1. 


1.5 Two spins | 

For two spins we consider a four dimensional Hilbert space which is a tensor 
product of the two dimensional Hilbert spaces of the one spin system. As the 
energy is given by the eigenvalues of the Hamiltonian, energy conservation 
needs that the eigenvalues of the sum matrix is equal to the sum of the 
eigenvalues of each Hamiltonian matrix. This property is obtained by the 
Kronecker sum ©y of two matrices which is defined by: 

A@ K B = A® I n + I m ® B, (52) 
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with 8) the tensor product with identity matrices / of dimensions n the first 
dimension of matrix B and m the second dimension of matrix A. 

Therefore the Hamiltonian of the two spin system is: 

= (k[ - fa • Si) ©A- (k ' 2 - fa ■ B 2 ) (53) 

This description is the most general possible for two independent spins in a 
B held. The description we have obtained for the wave function of two spins 
is easily generalized and the Hamiltonian of N spins is just the Kronecker 
sum of the one spin Hamiltonian: 

N 

^ l- (54) 

2—1 

This property shows the importance of the Hamiltonian of one spin which 
determines the dynamic of the whole free spin system. 


2 Measurements 


In NMR, the measurement of the spin system is obtained by the current 
induced in a tuned LCR circuit which corresponds to the photon emission 
when a spin switch from the excited state to ground state. It is well known 
that after a RF pulse an electrical current is induced in the LCR circuit, 
but even in absence of RF pulse, the presence of the sample can be detected 
by the nuclear-spin noise 0 . 0 . d, 0 . This signal has even been used for 
imaging sample without any RF pulse (TJ but in this case we will see that the 
presence of magnetic held gradients explains an increase of the NMR signal 
which permits to obtain images with a relatively good signal to noise ratio. 
In order to discuss the difference between these two cases, we calculate the 
probability of a spin to be in each state from equations (51a) and(51b) : 


x 2 x 2 = r ‘2 + sin ft 
XiX\ = r\ + sin ft 


rlA + sin ft + 2f±-r 1 r 2 B 
r\A + r 2 2 ^ sin ft - 2 ^r x r 2 B 


with 

A = sin ft($ - 1) 

B = (cos ft sin (0i - 02) + jy sin ft cos(0i - 0 2 ) 


(55) 
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which shows that the probability of a spin to be in one state is periodic 
around the value at t=0 and the period is Therefore the probabil¬ 
ity for a spin to switch from a state to the other depends on the value of 
A = + ujy + K 2 + K'jj X - For a free spin in a very homogeneous 

magnetic held A = K (remember that ui z = ujq — ca(0) = 0) and the proba¬ 
bility for a spin to change from one state to the other at a rate greater than 
the spontaneous oscillation between the two states is A j n presence of a 
RF magnetic held, we are in the case where cm = ujq and at resonance the 
probability to change from one state to the other is 1 because ojx, wy and 
K are small compared to cn 0 - 

2.1 Nuclear spin noise experiments 

After the hrst observation of random noise emission from nuclear spins by 
Sleator, Hahn, Hilbert and Clarke 0 , 0 using a superconducting SQUID, 
McCoy and Ernst [SJ have observed nuclear spin noise at room temperature 
using a conventional (300 MHz) NMR spectrometer. They compared the 
proton resonance of the methyl group in ethanol without any RF pulse just 
by co-adding 8000 spectra with a conventional one pulse spectra. The most 
surprising features of the spectrum without any RF is the negative devia¬ 
tion compared to the spectrum obtained by a 7t/2 pulse excitation and the 
amplitude of the signal which is about 10 s time smaller. Our elementary 
derivation retrieves this main characteristics. 

Without RF and no held inhomogeneity, the probability for a spin to emit 
a photon is given by —K/uj 0i as Ii = where K' is the kinetic energy 
of a proton K' = A- therefore the probability of a NMR signal detection is 

- — = 2 1CD 16 for protons in a 300 MHz magnetic held as in the McCoy 

and Ernst experiment. This probability is very small therefore in case of a 
great number of spins the statistical fluctuations around the mean values of 
r\ and r 2 gives a squared probability which is \/l 10 ” 16 ~ 1CT 8 . The minus 
sign gives the deviation from the probability to End the spin in a state at 
t — 0. In case of a RF magnetic held at resonance the probability is one and 
the deviation from state at t=0 is positive. The elementary description of 
the spin motion retrieves the main characteristics of a spin noise spectrum 
compared to the spectrum obtained with a 7 t/ 2 pulse. 
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T 7 =* • * ’" T - T ■ ■— 

Figure 1: Conventional gradient echo imaging sequence. 

2.2 NMR imaging experiments 

In imaging experiments, magnetic gradient fields are applied in all space 
directions as in the work of Muller and Jerschow [I] or in the three main 
direction as in phase encoding imaging thechnics. In this case a first gradient 
is applied in the X-direction with a RF pulse in order to select a slice with 
its thickness, after that a gradient is applied in a orthogonal direction in 
order give a phase coherence and after that a gradient is applied in the third 
orthogonal direction during signal acquisition. Figure [T| gives the schematic 
of one of the simplest imaging sequence. The first line corresponds to the RF 
channel where “om” is the synchronization of the lock-in detector with RF, the 
gaussien pulse shape is the red line and (o) corresponds to the acquisition of 
the current induced in the NMR receiver. The second line corresponds to the 
applied magnetic field gradients in blue for the Z (or slice) direction, in red for 
the Y (or phase) direction and in green for the X( or read) direction. When 
applying gradients the exact synchronization on the reference frequency ojq 
is mandatory to have the t = 0 values and the same X-direction for each 
phase encoding step or gradient orientation. Image A of figure [2] is the 
conventional NMR image of a test object composed of two Lego pieces in 
water with lg/L C 11 SO 4 and Agar/Agar 10 g/L. Images of 128x128 pixels 
was obtained on a Bruker 7T BS70/30 small animal imaging system. Image 
A of figure (2) was obtained with field of view of 58 mm corresponding to 
a magnetic field gradient of 2 10 -2 T/m, a slice thickness of 1.16 mm, and 
an echo-time of 3.9 ms with a repetition-time of 8 ms corresponding to an 
acquisition-time of Is 23ms. The RF pulse used was a Gaussian pulse of 15 
mW power corresponding to a flip angle of 5°. The signal to noise ratio in 


14 











Figure 2: NMR. images of a test object (Lego pieces): A with RF and B 
without RF pulses. 


image A is 18.1. Image B of figure [2] was acquired with the same parameters 
except that there was no RF pulse (Power excitation pulse <10 _9 PF), in this 
case the signal to noise ratio is 5.0. We see from equations (55) that NMR 
images can be obtained without RF magnetic held due to the fact that the 
probability to emit a photon by the spin system is not very different with or 
without RF because A without RF is A = y/u >| + u\ + Uy + K 2 + 2Kuj x 
and with RF A = \J oo\ + c o\ + Wy + K 2 + 2K(u>x + oonpcosoot ) + oo 2 RF . In 
effect in presence of magnetic held gradients A is not very different in both 
cases because oo RF is always small compared to oox, and the variation of 
ooz- An adequate discussion of the enhancement of the NMR signal with and 
without a RF pulse needs a precise description of the magnetic properties of 
the system under investigation and the information on spatial distribution 
of the spins, but the important result is that NMR imaging with good signal 
to noise ratio can be obtained without any RF magnetic held. Such imaging 
sequence can be of great medical interest for decreasing RF power deposition 
in a patient. 


Conclusion 

This work has shown that an elementary approach of the quantum descrip¬ 
tion of an ensemble of independent spin | can retrieve the magnetic resonance 
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basic experimental facts. This derivation completes the density matrix for¬ 
malism which is well adapted for describing all NMR experiments where spin 
coupling has important effects. Experiments using nuclear spin noise are dif¬ 
ficult to be explained by the matrix density formalism without introducing 
an extra probability factor and a stochastic operator as in the paper of Field 
and Bain |9J. A precise description of NMR quantum effects should partici¬ 
pate to the quantum mechanic fundamental aspects used in the development 
of quantum computing where NMR is still the only example of an experi¬ 
mental quantum computer working with more than ten qubits. The main 
result of our approach is to predict that the quantum interference effects are 
best detected in very homogeneous magnetic held without any RF pulse. 
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